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ABSTRACT 

Aims. We investigate the electromagnetic interaction of a relativistic stellar wind with small bodies in orbit around the star. 

Methods. Based on our work on the theory of Alfven wings to relativistic winds presented in a companion paper, we estimate the 
force exerted by the associated current system on orbiting bodies and evaluate the resulting orbital drift. 

Results. This Alfvenic structure is found to have no significant influence on planets or smaller bodies orbiting a millisecond pulsar. 
On the timescale of millions of years, it can however affect the orbit of bodies with a diameter of 100 kilometres around standard pulsars with 
a period P ~1 s and a magnetic field B ~ 10 8 T. Kilometer-sized bodies experience drastic orbital changes on a timescale of 10 4 years. 
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1. Introduction 

Accretion discs are expected to form at some phase of the evo- 
lution of neutron stars in a binary system, possibly giving birth 
to second generation planets. The interaction of the disc with 
the pulsar's wind constrains its extension and notably the value 
of its inner radius. The rotational power transferred to the disk 
by the magnetic field of a young neutron star (propeller effect) 
would rapidly disrupt a disc of plasma, put ting its inner radius 
beyond the light cylinder (lEksi et al.l l2005). Concerning dust, 
ablation by the pulsar's wind would extract protons and pro- 
vide them enough kinetic energy to c hase them awav (Ijonesl 
2007b . The Poynting-Robertson effect dBurns et al.ll 1979b may 
also act on the drif t of smaller (and isotherm al) particles such 
as grains and dust (ICordes & Shannonl2 008). This is consistent 
with infra-red observations of radio-pulsars showing that their 
inner radius is two or more orders of magnitude larger than the 
light cylinder radius (Jones 2008). 



On the other hand, small bodies, such as planetoids, aster- 
oids or comets may also orbit pulsars and it is expected that 
some of them occasionally f all be low the light cylinder. For 



Cord es & Shannon! d2008l) have shown that neutral, 



instance 

circumpulsar debris that enters the magnetospheres of neutron 
stars could disrupt current flows and account for some of the 
intermittency seen in radio pulsars. The neutral material of size 
~ 1 m or more can move toward the star under the influence of 



collisions and of the Yarkovsky effect (a net force induced by 
a difference of temperature between the illuminated and warm 
afternoon side and the night-time face). 

In the present paper, we study the influence of a pul- 
sar's wind on the trajectory of hypothetical larger bodies (1- 
100 km) through the action of the Poynting flux. It is based 
on Alfven wings, whose th eory is described in a joint paper 
dMottez & Hevvaertsfcoi lb . hereafter (M H1). The concep t of 
Alfven wings was initially developed by iNeubauei ( 1980b to 
explain the interaction of Jupiter with its satellite Io. It says that 
a conducting solid body embedded in a plasma flow slower than 
Alfven waves supports a system of electric currents carried by 
a stationary Alfvenic structure. When the wind is relativistic 
albeit slower than the total Alfven speed, the amplitude of this 
current can be estimated. In the present study, we estimate the 
force exerted by the associated current system on orbiting bod- 
ies and evaluate the resulting orbital drift. 



2. Alfven wings 

In this section, results demonstrated in the companion paper 
(MH1) and useful for the present study are recalled. 

The star is assumed to be spinning with an angular fre- 
quency Q„. For simplicity, the theory is developed for a pulsar 
with a magnetic dipole axis aligned with the rotation axis. The 
wind flow is assumed to be radial (r direction in spherical coor- 
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dinates) and the magnetic field in the companion's environment 
to be mainly azimuthal (0 direction). The wind is characterized 
by two invariants along its flow : the neutron star's magnetic 
flux *P, and the mass flux /, 



/ 



WW 
r 2 BL 



(1) 

(2) 



where vo is the unperturbed wind's velocity, B r Q is the radial 
magnetic field, jq is the wind's Lorentz factor, and p' {) is the 
proper rest mass density in the unperturbed wind's frame. The 
azimuthal magnetic field is given by, 



m — n r 



•Q.r 



(3) 



the approximation in the right hand side term being relevant at 
large distance from the light cylinder. 

The engine of the Alfven wings is the convection electric 
field associated to the wind that appears in the reference frame 
of the star's companion. The wing currents are generated by a 
potential drop U along the body of radius Rp, 



U = 2R P E = 



(4) 



where Eq - Vq-Bq, directed perpendicularly to the wind flow 
and to the magnetic field, is the convection electric field in- 
duced by the motion of the wind into the magnetic field. This 
potential generates a system of currents that flow along the 
companion, then in space into the plasma, in a direction which 
depends on the wind's magnetic field and the wind's velocity. 
The conducti vity of the plas ma part of this circuit has been 
evaluated bv lNeuba uer ( 1980) for the non-relativistic motion of 
Io in the Jovian magnetospheric plasma. In the ultra-relativistic 
wind of a pulsar, it can be approximated very simply by 



1 

fi c' 



(5) 



(See Eq. (67-69) of (MH1) for more details.) 

Then, adopting a simplifi ed geometry, it is possible to esti- 
mate the total electric current . INeubauerl (1 1 9 8 ) gives useful ex- 
pressions for the total current / flowing along an Alfven wing. 
Writing R P for the body's radius, he gets: 



I = 4(E -Ei)Rp-L A =4 



- Ej Rp T.a 



(6) 



The electric field set along the body, is caused by its iono- 
sphere or surface internal resistance. The Joule dissipation is 
maximum when £, = Eq/2. In our estimations, we shall use 
Neubauer's values for /. 

3. A magnetic thrust 

The above theory, because of the involved symmetries, de- 
scribes mainly what happens in space, far enough from the 
body. At closer distances, the plasma suffers compressive mo- 
tions and compressive MHD waves certainly have a non- 
negligible influence on the system. These waves propagate 



quasi-isotropically. Their amplitude decrases as the inverse of 
the distance to the body and they contribute to deflect the wind 
around it. Nevetheless, without entering into these considera- 
tion, we can still make a few inferences based on the theory 
of the Alfven wi ngs, as presented in the previous section. As 
Neubauerl dl980h . we can assume that the current associated 
to the wing is closed in the vicinity of the body (see Fig. [TJ, 
through its surface or its ionosphere. We can estimate (roughly) 
what force the wind exerts on it. 

The two Alfven wings carry a current that, in the two 
branches flowing along the body (perpendicular to the plane of 
FigHK generates a force density j x B. The two current systems 
flowing on each side of the body exert this force density in the 
same radial direction. We may expect the body to orbit near the 
equatorial plane of the pulsar. In this plane, at such a distance, 
the magnetic field direction is almost azimuthal, being perpen- 
dicular to the wind flow velocity. The sign of B^ Q depends on 
whether the magnetic moment of the neutron star is parallel or 
antiparallel to the rotation axis. Nevertheless, the force density 
j x B always has the same direction as the wind velocity. 

At first order, considering Eq. (0), Eq = Q.^/r and the 
force is expressed explicitely as a function of the distance r 
from the pulsar to the body as 



F = 2R P IB' P = 8 



EAR 



p y 



(7) 



The power E j dissipated by Joule effect along the ionosphere 
or in the body is maximized when the i nternal loa d matches the 
external one, that is, still according to Neubauer(1980), when 
Ej = Eq/2. In that case the force is 



4 ^q2vj/2 
fioc 2 r 2 



(8) 



On the night side of the body, this force tends to wipe out the 
ionosphere (if there is one), but on the day side, on the con- 
trary, it pushes the ionosphere towards it. The dynamics of 
this system is probably quite complex, but we can retain that 
there is a force pushing the body and/or its atmosphere away 
from its star. Maybe the atmosphere has been completely wiped 
out, and the day side of the body is ionised by the flux of X 
rays coming from the neutron star. Then, the current may flow 
along the dayside of the body's crust, directly pushing it away. 
Actually, it remains that, when the field is azimuthal and the 
dissipation is maximal, the force is radial and proportional to 
r~ 2 (Eq. ([8])). Therefore, it acts the same way as the gravita- 
tional force, and cannot have a secular influence on the orbit. 
When Ej does not vary with the distance r of the body as Eq 
does, according to Eq. (0, a fraction of the force is not of a 
Keplerian nature. For an azimuthal field, this force is however 
still central and causes nothing more than a periastron preces- 
sion. Since we are mainly interested in the evolution of the 
semi-major axis and of the eccentricity, we don't consider this 
case any further in this paper. 

But the unperturbed magnetic field also has a small radial 
component and the force is therefore not exactly central. Our 
estimate that B^ ~ fi' Q»r/c shows that the small angle between 
the magnetic field and the azimuthal direction is: 



5 = c/0,r. 



(9) 
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In the case of PSR 1257+12, at 1 AU, 6 = 2 x 1(T 6 . In the 
general context of a vacuum dipole wave, or a pulsar wind, its 
sign does not vary, and its amplitude decreases gently with the 
distance. As the force density j x B is perpendicular to the mag- 
netic field, the force is not strictly radial when the field is not 
strictly azimuthal. In spite of the small value of this angle, this 
azimuthal force component acts constantly in the same direc- 
tion. Therefore, this force can work. 

The tangential component of the force always has the same 
direction as the rotation of the neutron star. Therefore, if the 
planet's orbital angular momentum and the star's rotational 
spin are parallel (in the same direction), the j x B force con- 
tributes to its acceleration, and therefore it increases its semi- 
major axis and its eccentricity. The force modulus increases at 
smaller distance and the angle 5 also becomes larger. These two 
effects cause the tangential force to become stronger at closer 
distances from the star. 

From Eqs. (© and ([9j, a rough estimate of the tangential 
force is 

4 Rl^Cl* 



F, = F5 



f4>c 



(10) 



Let v or b ~ (GM./r) 1 ' 2 be the orbital velocity. The power asso- 
ciated to the work of F, is 



W, = F, Vo 



AG 111 Ml /2 R 2 J¥ 2 Q, 



r 7/2 
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Four planets have been discovered around two pulsars. 
Three of them orbit PSR 1277+12, with periods of the order 
of a few weeks. Their mass is comparable to the mass of the 
Earth (more details are given in (MH1)). Planets around pul- 
sars are expected to have been captured by the neutron star or 
to have orbited its progenitor before the supernova explosion. 
A capture would provide an initially large orbital eccentricity. 
Similarly, a body surviving a supernova explosion should be 
left after the event with a large eccentricity. However, the ec- 
centricities of the orbits of planets measured around pulsars are 
very small (see Table (2) in (MH1)). Could the tangential com- 
ponent of the j xB force associated to the planets' Alfven wings 
be an explanation ? 

Let us notice that a captured planet orbiting in the opposite 
direction to the star's rotation experiences a tangential compo- 
nent of the j x B force that tends to slow it down and/or to 
reduce its eccentricity. 

This force may also be exerted on smaller bodies, such as 
comets or asteroids. Here again, it is interesting to know how 
the orbits of such bodies would be influenced by their Alfven 
wings. This may be of importance for second generation plan- 
ets, which form (or not) after the supernova explosion from 
solid debris in the fall back accretion disc. 

4. Influence of the magnetic thrust on the orbit 

Let us now write down the equations of motion of an iso- 
lated body, orbiting a neutron star, under the action of the 
(Newtonian) gravitational force and the magnetic thrust. The 
magnetic thrust is decomposed into its radial component (Eq. 
(0) and orthoradial component. We roughly assume Ej to be a 



fraction of Eq. The orthoradial component is given by Eq. (flOt . 
The acceleration then is: 



M P 



% and 



M P 



D 

73' 



where C and D are constant factors, 



C = 



D 



4-R 2 p £l 2x ¥ 2 
fi c 2 M P ' 

4/? 2 Q„ > P 2 

\XQcMp 



(12) 

(13) 
(14) 



The radial force does not modify the body'orbit, which remains 
Keplerian, the star's mass M, being replaced by the slightly 
lower mass M, 

M - M, -— . (15) 
G 

The equations of motion then become 



GM 



d(r 



dt 



D 

~2 



(16) 
(17) 



These equations are those of a Keplerian motion, with a small 
correction induced by the right-hand side term of Eq. (TPTl . 
Therefore, we can consider that this motion is Keplerian in first 
approximation, and that the orbital elements a and e evolve 
very slowly. They are quasi-constant over an orbit. It is pos- 
sible to compute their slow average variations over an orbital 
period, < da/dt > and < de/dt >. One first needs to compute 
the instantaneous drifts da /dip and de/dcp, and the variations 
Aa and Ae over an orbit through an integration over <p from 
to 2n, considering the values of a and e to be constant in this 
sum and the motion to be purely Keplerian. The average values 
< da/dt > and < de/dt > are the variations Aa et Ae divided 
by the orbital period. The orbit of a Keplerian motion is repre- 
sented by the equation 



a(l -e 2 ) 
1 + e cos <b 



(18) 



The correction (proportional to D) is equivalent to a force that 
is tangential to the orbit. The in duced variations of a and e in 
this rather standard problem are (IMilani et al 1 ll987l) : 



da 2D (1 + e cos <p) 2 

~d(p ~ ~GM (\-e 2 ) 2 

de D e(\ + cos 2 ^) + 2cos< 

~chj> ~ GMa (1 -e 2 ) 

In one orbit, a changes by 

-In 



Aa 



r^da 
Jo dcp 



dip - 



AnD 1 + e 2 /2 



GM (I- e 2 ) 2 
The change of the eccentricity over one orbit is 



Ae 



r 2 " de 
Jo d(f> 



d<p — 



3nD 



GMa (1 -e 2 Y 



(19) 
(20) 

(21) 

(22) 
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The average changes of a and e over an orbital period are 



da 
dt 
de 
dt 



> = 



Aa 

— = 2a 
P 



D 



Ae_ 
~P~ 



y/GMa 5 
D 



2 + e 2 



2(1 - e 2 ) 2 



VGMa 5 (1 -e 2 ) 



(23) 



(24) 



These variation rates both have the same sign as D. We can 
see that for a prograde orbit, D > 0, a and e increase, the or- 
bit becoming more eccentric and distant. Therefore, the Alfven 
wing thrust tends to chase the body away from the star. A ret- 
rograd orbit evolves toward a circular shape with a decreasing 
semi-major axis. 

We now present a few numerical applications of Eqs. (1231 
and (l24l i. The basic numerical data about the pulsars and their 
companions can be found in Tables 1 and 2 of the joint pa- 
per (MH1). In table Q] of the present paper, we have written 
the corresponding yearly variations of a and e for the four 
know pulsar's planets, and for hypothetical pulsar's compan- 
ions. For large planets, the influence of the Alfven wing on the 
orbit is negligible; the Alfven wing could not explain why the 
orbits of the planets around PSR 1257+12 and PSR 1620-26 
have a small eccentricity. The orbit of small bodies (1 and 100 
km) orbiting a 10ms pulsar is also not significantly altered by 
the Alfven wing. On the contrary, the trajectory of a 100km 
body orbiting a Is pulsar, because of a much larger ambient 
magnetic field, can be significantly modified on a time scale 
of millions of years. The effect is significant in only 10 000 
years for a 1km sized asteroid. The orbits of asteroids orbit- 
ing in the sense of the neutron star's rotation spin increase in 
size and become more and more eccentric. Objects in counter- 
rotation (anti-parallel orbital angular momentum and rotation 
spin) would be quickly precipitated onto the neutron star. We 
can then expect that a young pulsar, even isolated (with no ac- 
cretion disc), can stimulate the in-fall of small bodies. 

In order to get a more precise idea about the Eqs. (1231241 . 
we have solved them numerically, using a fourth order Runge- 
Kutta algorithm, for a small set of initial orbits. We have 
solved the equation over a time span small enough to con- 
sider that the pulsar parameters that determine the parameter 
D are constant. Figure |2]shows the evolution of the semi-major 
axis of a 1 km size body, initially at a distance of 0.16 AU, 
for various values of the initial eccentricity. We have chosen 
D/ y/GM = +10 14 m 5 / 2 ■ s _1 that corresponds to the value given 
in Table Q] The sign + means that the body orbits in the sense 
of the pulsar' spin. We can see that, in accordance to what was 
said above, the semi-major axis increases significantly in a time 
scale of tens of thousand years. The increase is slightly larger 
for an initially large eccentricity. Figure |3]shows the evolution 
of the corresponding eccentricities. There is no variation for an 
initially null eccentricity, and the larger the initial eccentricity, 
the larger its further increase. 

Figure |4]shows the evolution of the semi-major axis for the 
same body in the case of counter-rotation. In less than 6000 
years, the body falls onto the star, or at least beyond the in- 
ner frontier between the wind and the magnetosphere (which 
is the limit of validity of the present calculations). In the in- 
ner magnetosphere, the magnetic thrust still acts onto the body 
(although differently), and a rapid fall onto the star is expected. 



5. Conclusion 

A planet orbiting around a pulsar develops a system of 
Alfven wings, caused by its interaction with the sub-Alfvenic 
Poynting-flux-dominated pulsar wind. A system of strong elec- 
tric currents is set, which exerts an ortho-radial force upon the 
planet that can, if the magnetic to mechanical energy coupling 
is efficient, have an incidence on the planetary orbit. 

The data in Table Q] show that the orbital drift of a planet 
around a millisecond pulsar is negligible. In particular, this ef- 
fect cannot be involved in an explanation for the very low ec- 
centricity of the four known planets orbiting a pulsar. But, for 
bodies with a diameter of a few kilometres orbiting around a 
P = 1 second pulsar, the drifts occur on time scales that are 
short in comparison with the time of evolution of an isolated 
pulsar. For bodies orbiting in the same direction as the star'spin 
(the orbital angular momenta being in the same direction as the 
star's angular momentum) this force tends to increase the semi- 
major axis and the eccentricity. For bodies in counter-rotation 
(the two angular momentum having opposite directions), this 
force tends to decrease the eccentricity and the semi-major 
axis, favouring the precipitation of the body onto the neutron 
star. If the bodies falling in t o the p ulsar's magnetosphere stud- 
ied by ICordes & Shannon! ([2008) fall under the influence of 
the Alfven wings, then these object were, before their fall, 
in counter-rotation with the the neutron star. Rocks and aster- 
oids in counter-rotation may come directly from the fall back 
after the supernova explosion. The y may as well come from 
recent collisions. As suggested by ICordes & Shannon! (|2008) 
the smaller bodies (about 1 metre) falling toward the star, and 
evaporating in its vicinity, may cause a momentary interruption 
of the pulsar's radio emissions. It is also possible that the fall 
of a larger object on the star's crust powers transient high en- 
ergy emissions, such as those observed with soft gamma-ray re- 
peaters. This might be particularly relevant for those assoc iated 
to pulsars with a "stand ard" magnetic fi eld (Rea et al.ll2010l) or 
with the Crab nebula dAbdo et alJuOl 11) . These aspects of the 
question remain to be explored 

More generally, the force caused by Alfven wings may have 
consequences on the dynamics of fall-back accretion discs, and 
onto the formation of second generation planets around pulsars. 
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Table 1. Electric potential drop, total electric current associated to the Alfven wing. Electrical energy Ej max dissipated in the 
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Fig. 1. Schematic view of an unipolar inductor. The unper- 
turbed wind's magnetic field Bo and velocity vo are almost, but 
not exactly, perpendicular. The electric field Eo created by the 
unipolar inductor is perpendicular to these two vectors; it in- 
duces an electric current (of density j) along the body. This 
current then goes into the interplanetary medium, forming two 
structures, each of them made of an outward and an inward 
flow. The current density j flowing along the planet is the cause 
of a j x B force density that is the topics of the present study. 
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Fig. 2. Evolution of the semi-major axis of a 1 km sized body as 
a function of time, under the influence of the magnetic thrust. 
The body orbits in the same direction as the pulsar' spin. The 
four curves are given for four different values of the initial ec- 
centricity. 



time (yr) 

Fig. 4. Evolution of the semi-major axis of a 1 km sized body as 
a function of time, under the influence of the magnetic thrust. 
The body orbits in the opposite direction to the pulsar' spin. 
The four curves are given for four different values of the initial 
eccentricity. 
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Fig. 3. Evolution of the eccentricity in the same cases as in Fig. 
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